Bulk and surface biaxiality in nematic liquid crystals 
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Nematic liquid crystals possess three different phases: isotropic, uniaxial, and biaxial. The ground 
state of most nematics is either isotropic or uniaxial, depending on the external temperature. Nev- 
ertheless, biaxial domains have been frequently identified, especially close to defects or external 
surfaces. In this paper we show that any spatially-varying director pattern may be a source of bi- 
axiality. We prove that biaxiality arises naturally whenever the symmetric tensor S = (Vn)(Vn) T 
possesses two distinct nonzero eigenvalues. The eigenvalue difference may be used as a measure 
of the expected biaxiality. Furthermore, the corresponding eigenvectors indicate the directions in 
which the order tensor Q is induced to break the uniaxial symmetry about the director n. We 
apply our general considerations to some examples. In particular we show that, when we enforce 
homeotropic anchoring on a curved surface, the order tensor become biaxial along the principal 
directions of the surface. The effect is triggered by the difference in surface principal curvatures. 
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PACS numbers: 

61.30.Gd - Oricntational order of liquid crystals 
61.30.Hn - Surface phenomena: alignment, anchoring, 
61.30. Jf - Defects in liquid crystals 



Nematic liquid crystals are aggregates of rod-like 
molecules. Early theories Q, y| used a single order 
parameter, the director, a unit vector pointing along the 
average microscopic molecular orientation. Most nematic 
phenomena fit well within the classical description. How- 
ever, the transition from ordered to disordered states es- 
capes the director theory. The classical microscopic de- 
scription of defects and surface phenomena yields unde- 
sired results as well. The order-tensor theory put forward 
by de Gennes 0, focuses on the orientational proba- 
bility distribution, and introduces the measures of the 
degree of orientation and biaxiality. Within this theory, 
a nematic liquid crystal possesses three different phases, 
which can be identified through their optical properties, 
since its Fresnel ellipsoid is closely related to the order 
tensor itself 6J. A isotropic liquid crystal is characterized 
by an isotropic order tensor, and optically behaves as an 
ordinary fluid. A uniaxial nematic possesses a unique op- 
tic axis. Its order tensor has two coincident eigenvalues. 
Finally, in a biaxial nematic the eigenvalues of the order 
tensor are all different, and the Fresnel ellipsoid possesses 
two optic axes. 

Within the Landau-de Gennes theory, the ground state 
may be either isotropic or uniaxial, depending on the ex- 
ternal temperature. However, biaxial domains have been 
predicted and observed, especially close to defects and ex- 
ternal boundaries. Schopohl and Sluckin yj analyzed in 
detail the biaxial core of a +i nematic disclination. More 
recent studied show that a biaxial cloud surrounds most 
nematic defects ||, and both analytic [UEJ and numeric 
[Til Il2j | asymptotical descriptions of biaxial defect cores 
have been derived. Other examples of defect-induced bi- 
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surface-induced ordering, 



axiality involve integer-charged disclinations 0, Q, 
and cylindrical inclusions [0|. The onset of surface bi- 
axiality is closely related to the presence of a symmetry- 
breaking special direction, which coincides with the sur- 
face normal 01 • Indeed, biaxiality has been predicted 
close to both external boundaries 0, 0] and internal 
isotropic-nematic interfaces pfl l2lj . 



In this paper we show that biaxiality effects are closely 
related to, but not exclusively confined to, the exam- 
ples above. In fact, within any spatially- varying director 
distribution, the director gradient itself breaks uniaxial 
symmetry about the director. We analyze in detail the 
structure of the elastic free energy density and come up to 
the result that, given the director distribution, it is pos- 
sible to predict the onset of biaxiality, to determine the 
direction of the secondary optic axis and to estimate the 
intensity of biaxiality effects. We then apply our general 
considerations to some specific examples, both within the 
bulk and close to an external boundary. We remark that 
we are not dealing with intrinsically biaxial nematic liq- 
uid crystals, that is systems in which the ground state 
itself becomes biaxial. Such systems, first observed by 
Yu and Saupe [2^, deserve a different treatment p3.l24[. 
since in them uniaxial symmetry is broken already at a 
molecular level. 



The paper is organized as follows. In Sections [I] and 
[n]we quickly review the order-tensor theory and the free 
energy density we aim at minimizing. In Section IlIII we 
derive and describe our main result, predicting a possible 
onset of biaxiality whenever the director is not uniform. 
In the following Sections IIVI and Ivl we apply the preced- 
ing results to some specific examples. In Section IVII wc 
collect and discuss our main results, while some Appen- 
dices collect the technical details of the proofs. 
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I. ORDER TENSOR 

The orientation of a single nematic molecule may be 
represented by a unit vector n£§ 2 , where S 2 is the unit 
sphere. Microscopic disorder is taken into account by in- 
troducing a probability measure f x : S 2 -» R + , such that 
/x(m) describes the probability that a molecule placed 
in x is oriented along m. The probability measure f x is 
even, since opposite orientations are physically equiva- 
lent. 

Nematic optics is determined by the variance tensor 
M = ( m ® m } , where the tensor product is defined 
in (| A8|) and the brackets denote averaging with respect 
to f x . By definition, M is symmetric and semidcfinitc 
positive. Since, in addition, the trace of M is equal to 1, 
we define the traceless order tensor Q = M — 1 1, where 
I is the identity. 

We label the nematic as isotropic when all the eigen- 
values of Q coincide, which implies Qi SO = 0. When 
at least two eigenvalues are equal, the nematic is called 
uniaxial. Simple algebraic manipulations allow to write 

Qu„i = s (n ® n - i ij , (1) 

The scalar parameter s is the degree of orientation |25j |. 
while the unit vector n is the director. The eigenvalues 
of Quni ar e |s (associated with n) and — |s (with mul- 
tiplicity 2). The director is then the eigenvector associ- 
ated with the different eigenvalue. Equivalently, n could 
be also identified as the eigenvector associated with the 
eigenvalue whose sign is different from the other two. 

When the eigenvalues of the order tensor are all differ- 
ent, the nematic is labeled as biaxial. In this general case, 
we can use the above remark, and still identify the direc- 
tor as the eigenvector of Q whose eigenvalue has a dif- 
ferent sign with respect to the other two. This definition 
may induce an artificial director discontinuity whenever 
the intermediate eigenvalue crosses 0. In turn, it yields 
an operative definition that works well when the order 
tensor is possibly biaxial, but however close to being uni- 
axial. Once we have introduced the director, we again 
define the degree of orientation s = |A n , where A n is the 
eigenvalue associated with n. The other two eigenval- 
ues A± can be finally written in terms of the degree of 
biaxiality (3: A± = — ^ s ± (3. As a result we obtain 

Qbia = s ® n - i ij + (3 (e+ ® e+ - e <g> e_) . (2) 

The sign of (3 is unessential, since it only involves an 
exchange between e + and e_. The degree of biaxiality 
does always satisfy \(3\ < ^\s\. Indeed, when \j3\ = ||s| 
one of the eigenvalues vanishes, and greater biaxiality 
values would in fact announce an abrupt change in the 
director (and in the degree of orientation as well) . 



II. FREE ENERGY FUNCTIONAL 

Equilibrium states of nematic liquid crystals are iden- 
tified as extremals of the free-energy functional whose 
density, in the absence of external fields, comprises two 
terms 

*(Q,VQ) = * cl (Q,VQ) + vI/ LdG (Q). (3) 

Though all the calculations we report could be re- 
peated in a more general framework, we will adopt the 
1-constant approximation for the elastic contribution ty c \ 

*d(Q,VQ) = y|VQ| 2 , (4) 

where K is an average elastic constant. 

The Landau-de Gennes potential 'I'LdG is a tempera- 
ture-dependent thermodynamic contribution that takes 
into account the material tendency to spontaneously ar- 
range in ordered or disordered states: 

*LdG (Q) = A tr Q 2 - B tr Q 3 + C tr Q 4 . (5) 

The material parameter C must be positive to keep the 
free-energy functional bounded from below. The poten- 
tial (J5J depends only on the eigenvalues of Q, and penal- 
izes biaxial states j2Q. Insertion of into (J5J returns 

*Ld G (s, P) = ~ (Cs 4 - Bs 3 + 3 As 2 ) 

+ I (6Cs 2 + 9Bs + 9A) 1 + 2C/3 4 . (6) 
9 

Let a — 3A/(Csq). The absolute minimum of \&LdG is 
located at the uniaxial configuration (sq > 0,/3 = 0), 
provided 

ae[-2,l] and B = § Cs (a + 2) . (7) 

When looking for minimizers of the free energy func- 
tional, we take into account that Landau-de Gennes' con- 
tribution usually dominates the elastic one. This approx- 
imation holds as long as we do not get too close to a ne- 
matic defect. Indeed, experimental observations confirm 
that neither s nor (3 depart easily from their preferred 
values (so, 0). 

We then envisage a two-step minimization. In the first 
step (s,/3) are constrained to their optimal values. Min- 
imization proceeds exactly as in Frank's director theory 
and yields an optimal distribution n(r). In the second 
step, we fix the director distribution and determine the 
perturbative corrections it induces in the optimal values 
of the scalar order parameters. As a result, we prove 
that non-uniform director configurations may induce a 
nonzero degree of biaxiality, and a reduction in the de- 
gree of orientation. As a by-product we determine how a 
non-zero director gradient breaks the local axial symme- 
try induced by the director, and which direction is chosen 
by most molecules (among those orthogonal to n). 
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III. BULK BIAXIALITY 

We collect in the present section our main result. In 
order to ease the reader we defer most of the technical 
proofs to the appendices below. We assume that a spe- 
cific director distribution n(r) has been determined by 
minimizing Frank's free-energy functional, constrained 
by suitable boundary conditions. The director distribu- 
tion may also take into account the effects of any possible 
external field. 

In Appendix lAl we prove the following decomposition 
for the director gradient 

Vn = A2 e 2 £S> e 2 + A3 e 3 g> e 3 + ( curl nAnjgn 
+ \ (n-curln)W(n), (8) 

where W(n) denotes the skew tensor associated with n 
(see Appendix 0). Furthermore, {A 2 ,A 3 }, {e 2 ,e 3 } are 
respectively the eigenvalues and eigenvectors of the sym- 
metric part of G = Vn— (Vn)n(g)n, the third eigenvector 
of sym G being n, with null eigenvalue. We remark that 



the term within brackets is maximized when {e + , e_} co- 
incide with the two eigenvectors of S that arc orthogonal 
to n. If we denote by the correspondent eigen- 

values, the linear term becomes simply proportional to 
— /!-). We thus arrive at the following result. Con- 
sider the symmetric tensor S = (Vn) (Vn) T . It always 
possesses a null eigenvalue (with eigenvector 11). When- 
ever its other two eigenvalues do not coincide, biaxial- 
ity is naturally induced in the system, and the optimal 
eigendirections of Q coincide with those of S. 

If we take into account expression (JSJ for Vn, we can 
give the eigenvalue difference (/i + — the following 
expression (see (lA7fl ) 

(ti + -^f = {cl-cl + >^-\lf 

+ 4 (c 2 c 3 + ±c„(A 3 - A 2 )) 2 , (11) 

where curln = c n n+c 2 e 2 +c 3 e 3 , and A 2 , A 3 are as in JSJl. 
In the following sections we will apply the above results 
to some practical situations, in order to better interpret 
their implications. 



div n = tr Vn = A 2 + A3 



(9) 



IV. SPLAY, BEND AND TWIST BIAXIALITY 



Let S be the symmetric tensor S = (Vn) (Vn) T . By 
virtue of IjAlj! the director n is an eigenvector of S (with 
null eigenvalue). In Appendix[5]we prove that the elastic 
free energy density may be given the following form 



*oi = K 



±|Vs| 2 + |V/3| 2 +s 2 |Vn| 2 



T e I 2 



+/3 2 (|Vn| 2 + 4|(Ve + ) 
-2s/?(e + • Se+ - e • Se 



(10) 



Let us analyze in detail the different terms appearing in 
(|10fl . The first two terms are trivial, since they sim- 
ply penalize spatial variations of the scalar order pa- 
rameters. They remind that, even in the presence of 
spatially-varying preferred values (s opt (r), /3 opt (r)) , the 
equilibrium distribution may not imitate the optimal val- 
ues. The third term is proportional to s 2 | Vn| 2 . This 
term has been already extensively studied |25t 12^ . Its 
net effect is a decrease in the degree of orientation in 
places where the director gradient is most rapidly vary- 
ing. In particular, it strongly pushes the system towards 
the isotropic state s — when the director gradient di- 
verges. The second-last term is proportional to f3 2 . Since 
it is positive-definite, it simply enhances the character of 
(3 = as optimal biaxiality value. Thus, were not for the 
final term we will next consider, biaxiality would never 
arise naturally in a nematic liquid crystal. 

The last term in (|10f> is linear in /3. It shifts the optimal 
biaxiality value away from (3 = 0. In order to minimize 
the complete free energy density it is worth to maximize 
the multiplying factor depending on S. This condition 
determines the directions {e + ,e_} in which the order 
tensor Q is pushed to break uniaxial symmetry. Indeed, 



A. Pure splay 

We begin by considering the splay field n(r) = e r , 
where e r is the radial unit vector in cylindrical co- 
ordinates. If we complete an orthonormal basis by in- 
troducing the tangential and axial unit vectors eg,e z , 
standard calculations allow to prove that the relevant 
fields are given by 



Vn = - eg 

r 



eg 



curl n = 



sym G = G = Vn = - e 



e 



1 



S = (Vn) 2 = — eg ®eg 



(Ve, 



= (Ve 8 ) r e 2 = 



(12) 



Thus, /!+ = r 2 , jit- 
density is given by 



K 



1 



|Vs 



= 0, and the elastic free energy 

l2 , (s-er 



|V/3|- 



(13) 



Biaxiality favours the tangential direction with respect 
to the axial direction. The r~ 2 factor implies that biax- 
iality (and the degree of orientation decrease as well) is 
expected to show close to the symmetry axis. Figures 3 
and 5 of 01 exactly confirm this result. 



B. Pure bend 

We again consider the same cylindrical coordinate 
frame above, and analyze the bend field n(r) = eg. We 
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obtain 



Vn = — e r ® eg 
r 

sym G = G = 



curln = - e 2 

r 



= (Ve r ) T e z = 



(14) 



Again, /x+ = r , /i_ = 0, and the elastic free energy 
density can given exactly the same expression <|13[1 . Bi- 
axiality now favours the radial direction, and again con- 
centrates close to the (disclination) symmetry axis. 



Pure twist 



We now introduce a Cartesian frame {e x ,e y ,e z } and 
consider the twist field n(r) = cos kz e x + sin kz e y . If we 



introduce the unit vector n j_ (r) 
we obtain 



v 

sin kz e,, + cos kz e 



Vn = k n± ® e 2 ==> curl n = — k n 
G = Vn => sym G = k sym nj_ <E> e 
S = k 2 n± (g> nj_ 



(Ve+) J e_ = (Vru) J e z =0 



(15) 



We now have //+ = k 2 , \i- = 0. Biaxiality favours nj_, 
that is, the (x, y) plane, with respect to the transverse 
direction e 2 . The elastic free energy density does again 
coincide with QliSfl. with only a k 2 replacing the r~ 2 fac- 
tor. However, this coincidence must not induce to guess 
that 'J'ei does always depend on s and /3 only through the 
combination (s — /3), as we will evidence below. 



D. Third-dimension escape 

All the model cases analyzed above share a peculiar 
property. Indeed, the fields n(r) considered are all pla- 
nar, in that they are all orthogonal to a fixed direc- 
tion (z-axis). Therefore it is of no surprise that we 
always find some degree of biaxiality which penalizes 
the avoided direction. We now consider a less triv- 
ial example: the escape in the third-dimension. This 
field was first determined by Cladis and Kleman [2^ 
as an everywhere continuous director field able to ful- 
fill homeotropic boundary conditions on a cylinder. Let 
n(r) = cos (j)(r) e r + sin</>(r) e z be the director field, and 
let nj_ (r) = — sin </>(r) e r + cos 4>{r) e z . We obtain 



Vn = — 

1 

curl n = 



eg 



cos m ea 



n 1 



r , r , cos 
sym (j = Lr — eg 



e e 



t>' sin <p nj_ (g> nj_ 



cos 2 <p 
a = — 7T- ea 



eg 



' nj_ (g) nj_ 



(16) 



Expression i|16[l for S shows that, within the order ten- 
sor Q, either n^ or eg may be preferred, depending on 
whether <fi' 2 is greater or smaller than cos 2 tfi/r 2 . This 
result turns out to be particularly challenging, if we con- 
sider that in Cladis-Kleman's escape in the third dimen- 
sion the tilt angle <p is given by 



. . 7T r 
(r) = 2 arctan — 



R 



(17) 



A simple calculation allows to show that i|17|) implies 
4>' 2 = cos 2 4>/r 2 . Thus, the third-dimension escape turns 
out to be one of the few spatially-varying director fields 
which do not induce any biaxiality. The elastic free- 
energy density in Cladis-Kleman's third-dimension es- 
cape is given by 



1 = K ( ~ iVa 



|V/3| : 



(r 2 + i? 2 ) 2 



4(i? 4 



r 2 (r 2 + R 2 ) 



V. SURFACE BIAXIALITY 



(18) 



In this section we estimate the degree of biaxiality in- 
duced by an external surface on which strong anchor- 
ing is enforced. We consider separately the cases of 
homeotropic and planar anchoring. Differential calcu- 
lus formulae that turn out to be useful for both cases are 
collected in Appendix IO 



A. Homeotropic anchoring 

We first assume that the surface director is parallel to 
the unit normal v to a given (smooth) surface E. We 
also assume that the director keeps its normal direction, 
at least in a thin surface slab. To be more precise, we pa- 
rameterize bulk points through a coordinate set (tt, 
such that 



P(u, v, £) = Pe(u, v) + £ u(u, v) 



(19) 



where Py. is the projection of P onto E, £ is the distance 
of P from the fixed surface, and v is the unit normal at 
Pe (see Figure QJ. If E if smooth, the coordinate set is 
well-defined in a finite neighborhood of E. We assume 
that n(P(u,v,g)) = n(P s (w,u)) = v(u,v). Then, Vn 
turns out to be closely related to the curvature tensor. It 
is symmetric and can be written as (see <|C2ll ') 



Vn = -- 



ei ® ei — 



l<2 



e 2 ® e 2 



(20) 



where {k\, K2} and {ei, e 2 } denote respectively the prin- 
cipal curvatures and principal directions at Ps. From 
them we obtain curln = 0, G = Vn, and 



■ ei ® ei 



(1 - 



■ e 2 ® e 2 



(21) 
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FIG. 1: Geometric setting for the surface parametrization 
introduced in the text. 



Equation (|21|l shows that biaxiality arises naturally close 
to an external surface where homeotropic anchoring is en- 
forced. This effect is triggered by the difference between 
the principal curvatures. More precisely, the tangent di- 
rection preferred by the order tensor is the one along 
which the surface curves more rapidly. Close to a sym- 
metric saddle, where ni — —K2, the denominator of l|21() 
induces biaxiality along the direction which is convex to- 
wards the side occupied by the liquid crystal. 



B. Planar anchoring 

When planar anchoring is enforced on a curved sur- 
face, it is natural to assume that the chosen direction 
coincides with one of the principal directions along E. 
We then keep the same notations as above, and assume 
n (P(u,v ,£)) = n(Ps(w, v)) — ei(u, v). When this is the 
case, by (IC2|) we obtain 







and 



u <g> ei 



curl n = 



(22) 



Thus, in the presence of planar anchoring, biaxiality 
arises whenever the curvature along the prescribed di- 
rection is different from zero. When this is the case, the 
biaxiality direction coincides with the unit normal. 



VI. DISCUSSION 

We have shown that any spatially- varying director dis- 
tribution may induce the onset in biaxial domains even 
in nematic liquid crystals whose ground state is strictly 
uniaxial. In particular, in Section II I II we have stressed 
the crucial role played by S = (Vn)(Vn) T . The tensor 
S, which is symmetric and positive-semidefinite by con- 
struction, possesses always a null eigenvalue, with eigen- 
vector n. Equation (|1L)|> shows that biaxiality arises nat- 
urally whenever the other two eigenvalues of S are dif- 
ferent. Then, equation shows that such eventuality 



is closely related to the vector curl n and the eigenvalues 
entering in the decomposition (JSJ of the director gradient. 

In Section llVl we have applied the considerations above 
to some model cases. As it could be easily predicted 
the pure splay, bend, and twist fields, being all pla- 
nar, exhibit some degree of biaxiality which privileges 
the director plane over the orthogonal director. A less 
trivial result is that there are spatially-varying direc- 
tor configurations that do not induce biaxiality at all. 
Cladis-Kleman's escape in the third dimension yields a 
unexpected example of this phenomenon. Section an- 
alyzes the onset of surface biaxiality both in the case of 
homeotropic and planar alignment. In the former case, 
biaxiality is ruled by the difference between the princi- 
pal curvatures along the surface. In the latter, only one 
curvature counts, and more precisely the one along the 
prescribed direction in the tangent plane. 

To conclude our analysis we want to give a numerical 
estimate of the magnitude of the biaxiality phenomena 
we are predicting. In all nontrivial cases, the free-energy 
density will contain a 0(/3)-term, which triggers the biax- 
iality onset. To obtain a rough estimate, we can neglect 
the 0(/3 4 )-term in * L dG, and the 0(/3 2 )-term in \& cl , both 
with respect to the dominant 0(/3 2 )-term, appearing in 
^LdG- When this is the case, the (local) preferred value 
of (3 may be obtained by minimizing the function 



g (j3) = - (QCs 2 + 9Bs + 9A) (3 2 - 2Ks(3(^ 



2(2 + a)Csll3 2 
2Ks 



£ 2 



(23) 



where we have replaced s 
matic coherence length 



£ 2 

Sn 



So and introduced the ne- 



K 



Cs (2 + a) ' 



(24) 



The (local) optimal value of the degree of biaxiality is 
then 



/3 op t « |fn (fi+ - M-) 



(25) 



Though /3 opt may vary from point to point, we have to 
keep in mind that in general the equilibrium configura- 
tion will not coincide with /3 opt because of the |V/3| 2 - 
term, and the boundary conditions. To make an explicit 
example, let us consider a nematic cylindric capillary of 
radius R, with homeotropic conditions enforced at the 
surface. Then, the difference between the eigenvalues of 
S at the surface is R~ 2 and the surface biaxiality is of 
the order of (£ n /R) 2 . Since the nematic coherence length 
hardly exceeds the tenths of a /xm, we obtain /3 pt < 10~ 2 
for a /im-capillary. The scenario changes completely close 
to a nematic defect, where at least one of the eigenvalues 
of S diverges. Both the pure-splay and the pure-bend 
examples above yield (fi + — /i_) = r~ 2 , which implies 
/3 op t ~ (£,n/ r ) 2 - The biaxiality cloud cannot be neglected 
if we come too close to the defect. 
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APPENDIX A: DIRECTOR GRADIENT 

In order to characterize the tensor Vn we begin by 
noticing that 

(Vn) T n= iv(n-n) =0 , (Al) 

since n is a unit vector. Thus, 

( Vn)n = ( sym Vn + skw Vn) n 

= i (Vn)n+ i curlnAn, (A2) 

and (Vn)n = curln A n. Let G = Vn — (Vn)n n. For 

any vector v, 

(skw G)v = (skw( Vn) - ± [( Vn)n n - n (Vn)nj) v 
= 5^ curln - n A (curln A n)J A v (A3) 
= |(n • cuiin)n A v = |(n ■ curln)W(n)v , 

where W(n) denotes the skew tensor associated with n, 
that is the tensor such that W(n)v = n A v for any v. 
Thus, 

Vn = sym G + \ (n • curl n) W(n) + ( curl n A n) n 
= A 2 e 2 e 2 + A 3 e 3 e 3 + |(n • curln)W(n) 
+ ( curl n A n) ® n , ( A4) 

where {A 2 , A 3 } and {e 2 , e 3 } are respectively the eigenval- 
ues and eigenvectors of sym G. The eigenvectors {e 2 , e 3 } 
are orthogonal to n, since <|A3ll implies 

(sym G)n = Gn — (skw G)n = (Vn — (Vn)n ® n)n 

- ~(n-curln)W(n)n = . (A5) 

Let us consider the symmetric tensor S = (Vn)(Vn) T . 
In view of the crucial role it plays in inducing biaxiality 
we now analyze it in more detail 

S = [A 2 e 2 e 2 + A 3 e 3 ® e 3 + |(n • curln)W(n) 

+ ( curl n A n) n] [A 2 e 2 ® e 2 + A 3 e 3 ® e 3 

— h(n ■ curln)W(n) + n ® (curln A n)] 

= A 2 e 2 ® e 2 + (A 2 - A 3 )(n • curln) sym(e 2 ® e 3 ) 

+ A 3 e 3 ® e 3 + |(n • curln) 2 (I - n <8> n) 

+ ( curl nAn) ® ( curl n A n) . (A6) 

Let {0,/i+,/i_} be the eigenvalues of S. The onset of 
biaxiality depends whether they latter two are equal or 
not. Let curln = c„ n + c 2 e 2 + c 3 e 3 . From (|A6|) we 
obtain 

+ 4 (c 2 c 3 + ±c n (A 3 - A 2 )) 2 . (A7) 

We finally remind that, given two vectors u,v, the ten- 
sor product (u ® v) is defined as the second order tensor 
such that 

(u ® v)a = (v • a) u for any vector a . (A8) 



APPENDIX B: ORDER TENSOR GRADIENT 

Let us differentiate equation J2J. We obtain 

VQ = (n ® n - ±l) ® Vs + s (Vn n + n Vn) 

+ (e+ ® e+ - e ® e_) ® V/3 + /5(Ve+ e+ 
+ e + Ve + - Ve_ e - e ® Ve_) , (Bl) 

where, given the second order tensor L and the vector u, 
(L u) is defined as the third order tensor such that 

(L u)a = La u for any vector a . (B2) 

When computing the square norm of VQ, we can make 
extensive use of the property (|A1|) and also take into 
account that 

uv = =^ (Vu) T v = -(Vv) T u . (B3) 
As a consequence, we obtain 

|VQ| 2 = | |Vs| 2 + 2 |V/3| 2 + 2s 2 |Vn| 2 

+ 2/3 2 (|Ve+| 2 + |Ve_| 2 + 2|(Ve+) T e_| 2 ) 

- 4s/?(|(Vn) T e + | 2 - |(Vn) T e_| 2 ) . (B4) 

We can further simplify expression (|B4|I if we consider 
that 

|Ve+| 2 + |Ve_| 2 = |(Ve + ) T | 2 + |(Ve_) T | 2 

= |(Ve + ) T n| 2 + |(Ve + ) T e_| 2 

+ |(Ve_) T n| 2 + |(Ve_fe + | 2 

= |Vn| 2 + 2|(Ve+) T e_| 2 (B5) 

and that 

(Vn) T u| 2 = (Vn) T u • (Vn) T u = u Su , (B6) 

provided we define S = (Vn)(Vn) T . By using {B5j-{B6j) 
it is immediate to give (|B4ll the expression quoted in HUfl . 

APPENDIX C: CURVATURE TENSOR 

Let E be the smooth surface, which bounds the sys- 
tem we are interested in. Let u be the unit normal, 
everywhere pointing in the direction of the bulk. We 
parameterize points in the bulk through a coordinate set 
(u, v,£) such that 

P(u,v,Z) = PE(u,v)+tv(u,v), (CI) 

where Py, is the projection of P onto E, and £ is the 
distance of P from the same surface. Such coordinate set 
is well-defined in a finite neighborhood of E. 

Let us consider the vector field everywhere defined as 
t(P(u,v,£)) — v(Pt:(u,v)). The second-order tensor 
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Vr is symmetric. It generalizes the curvature tensor 
V^i/, which is defined only on the tangent bundle of S. 
The eigenvectors of Vt coincide with those of the cur- 
vature tensor, and are thus the unit normal (with null 
eigenvalue) and the (tangent) principal directions on S. 
If we introduce {ki,K2}, the principal curvatures on S, 
and their corresponding eigenvectors {ei, 62}, we have 

Vr = -— — -ei (g> ei - — ^— - e 2 <g> e 2 (C2) 
1 - «i? 1 - K2C 

_ Ki , _ K2 

Vei = 1/ ® ej and Ve 2 = 1/ <g> e% . 
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